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We consider the form of the current-voltage curves generated when tunneling spectroscopy is 
used to measure the energies of individual electronic energy levels in nanometer-scale systems. We 
point out that the voltage positions of the tunneling resonances can undergo temperature-dependent 
shifts, leading to errors in spectroscopic measurements that are proportional to temperature. We 
do this by solving the set of rate equations that can be used to describe electron tunneling via 
discrete quantum states, for a number of cases important for comparison to experiments, including 
(1) when just one spin-degenerate level is accessible for transport, (2) when 2 spin-degenerate levels 
are accessible, with no variation in electron-electron interactions between eigenstates, and (3) when 
2 spin-degenerate levels are accessible, but with variations in electron-electron interactions. We also 
comment on the general case with an arbitrary number of accessible levels. In each case we analyze 
the voltage-positions, amplitudes, and widths of the current steps due to the quantum states. 
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I. INTRODUCTION 

Nanometer-scale single-electron tunneling transistors 
can now be fabricated in which electron flow occurs 
through a discrete spectrum of well-resolved quantum 
states. This has been achieved in devices incorporat- 
ing semiconducljiHg quantum dots, metal nanoparticlcs, 
and molecules .Mm In a transistor geometry, the source- 



drai n voltage V and the gate voltage V g Cctn be adjusted 



to achieve the simplest ca3c that electron flow occurs 

just through a single quantum state. As V and V g are 
changed, additional excited electronic states may also be- 
come energetically accessible for tunneling, providing al- 
ternative channels for current flow. In this regime, the 
tunneling processes can become quite complicated, due 
to the many combinations of non-equilibrium states that 
may be excited during tunneling, and the possibility of 
relaxation between these states. 

As long as the tunnel-barrier resistances are much 
greater than h/e 2 and internal relaxation is negligible, 
the currents traveling via any number of energetically- 
accessible states can be analyzed in a sequential- 
tunneling picture using a rate-equation approach. The 
general procedure for completing this type of analysis 
has been outlined previously, for example in Ref. |],|]]|. 
Our purpose in this paper is to present the solutions of 
this model for selected simple cases important for ana- 
lyzing experiments on non-magnetic islands, and we de- 
scribe several previously-unappreciated consequences of 
the model that explain recent observations. Whenever 
more than a single (non-spin-degenerate) quantum state 
is accessible for tunneling, we show that the voltage- 
positions of the tunneling resonances can become tem- 
perature dependent (shifting proportional to T). For 
the important case of tunneling via one spin-degenerate 
quantum state, we derive the full form of the tunneling 
current as a function of V, V g , and T. This provides 
simple exact solutions for the voltage shift, resonance 
width, and current amplitude, thereby improving upon 



an approximate approach used previously. When multi- 
ple spin-degenerate states participate in tunneling, effects 
of non-equilibrium excitations and variations in electron- 
electron interactions can lead to additional shifts and 
broadening of the tunneling resonances. The computer 
code which we use for calculating the tunneling current 
in the general case with an arbitrary number of accessible 
quantum states is available electronically in both Math- 
cmatica and C formats at ittp : //www. ccmr . Cornell . 



edu/~ralph/ 



This paper is organized as follows: In section II we 
review the general procedure for calculating tunneling 
currents in the rate-equation approach. We discuss the 
physical assumptions under which this approach is accu- 
rate, and we explain our notation. In section III, we solve 
the simplest non-trivial case, in which current flow occurs 
by means of tunneling via a single spin-degenerate quan- 
tum level. In section IV we then extend this discussion 
to the case of tunneling via 2 or more spin-degenerate 
levels, and we describe several experimentally-relevant 
consequences of the rate-equation model for an arbitrary 
number of accessible states. In section V, we consider ef- 
fects of fluctuations in electron-electron interactions that 
can occur when current flow generates non-equilibrium 
electronic states, and we explain how these effects can 
produce additional shifts and can also broaden the mea- 
sured tunneling resonances. 



II. RATE-EQUATION CALCULATIONS OF 
CURRENT FLOW 



We are interested in calculating the tunneling cur- 
rent via a non-magnetic single-electron transistor in the 
regime where the discrete quantum states in the transis- 
tor island are well resolved. The circuit under consider- 
ation is shown in Fig. [II which illustrates the definitions 
of the bias voltage V and the gate voltage V g . We will 
limit our discussion to the conditions under which the 
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FIG. 1: Circuit schematic defining the bias voltage V, the 
gate voltage V g and the capacitances Ci, C r and C g . <f> and 
Q are the potential and the total charge of the island. 



energy levels are best resolved: (a) ksT is smaller than 
the level spacing, (b) the level spacing is much smaller 
than the Coulomb charging energy of the transistor island 
e 2 /(2Cs), where Cs is the total capacitance of the island, 
(c) the tunnel barriers have resistances 3> h/e 2 so that 
cotunncling processes may be neglected and the tunnel- 
ing current is accurately described by lowest-order per- 
turbation theory, and (d) UbT is larger than the intrinsic 
lifetime broadening of the quantum states. In parts of the 
discussion, in order to simplify the notation, we will also 
assume that electron interactions are sufficiently weak 
that many-body eigenstates |a) are well-approximated 
as single Slater-determinants specified by the occupation 
of a set of single-electron states i: \ce) = {rii}. Our 
primary goals are to study the effects on current flow 
of non-equilibrium electronic excitations and electron- 
electron interactions. Non-equilibrium excitations can be 
suppressed when the rate of internal relaxation within 
the transistor island is large compared to the tunnel- 
ing rate. However, measurements on metal nanoparti- 
cles indicate that the relaxation rate is generally com- 
parable ton or slower than the tunneling rate in realistic 
samples. uu Therefore we will generally neglect internal 
relaxation effects entirely, limiting ourselves to noting the 
ways in which internal relaxation will produce qualitative 
changes to the results. 



A. Energy of the eigenstates 

In general, the quantum-mechanical electronic states 
within the transistor island can be complicated correlated 
many-electron eigenstates. The energy of any state can 
be written as a sum of three terms: 



E = Er 



E 



K 



(1) 



the terms being respectively the electrostatic or 
"Coulomb" energy, the kinetic energy, and the fluctu- 
ations in the electron-electron interactions. Notice that 
the mean-field contribution of the electron-electron in- 
teractions is the same as the electrostatic energy Eq- 
Therefore Ej accounts only for the level-to-level fluctu- 
ations in these interactions. 



1. Electrostatic energy 

The electrostatic energy will in general depend on the 
charge of the island as well as on the applied voltages V 
and V g . However, what matters for calculations of elec- 
tronic transport are energy differences as electrons make 
transitions between the transistor island and the leads. 
We can select our zero of energy (or, equivalently, the 
reference electrostatic potential) for convenience, and we 
will do so in a way that makes the energy of the eigen- 
states on the island independent of V and V g . The con- 
sequence is that the Fermi energies in the leads will shift 
with V and V g . To be specific, we choose the reference 
electrostatic potential such that J^k Ck4>k = 0, where Ck 
and 4>k are the capacitance of the island to the &-th lead 
and the electric potential of the /c-th lead, and the sum 
extends over the three leads. Using this reference, the 
charge Q in the island is related to its potential <f> by 
Q = Cy,4>- In calculating the energy required for a tun- 
neling transition, we must consider the work done. The 
tunneling of charges 8Qi and 8Q r from the island to the 
left and right leads requires a work 



SW = {4>i - 4>)8Q l + {4> r - 4>)5Q r 
= 4>i6Qi + 4> r 8Q r + —QSQ 



S <j>lQl + 4>r 



2CV 



(2a) 
(2b) 

(2c) 



where Qk is the total charge that has tunneled into lead 
k (Note SQ = -(SQi + 5QM From Eq. (|) it follows 
that the electrostatic energji!3 of the island is 



2C S 



(3) 



and the effective Fermi energies of the leads can be writ- 
ten as E£ = e<j>k , where e is the electron charge, including 
its sign. To be explicit, 



1C r + Cg 



V 



9 



2C i + C a . v _ e ^ Vc 



2C S 



(4a) 
(4b) 



Since the charge of the island varies only by multiples 
of e, we can write it as Q = Qq + ./Ve where Qq is a 
background charge and N the number of electrons in the 
island. The electrostatic energy is then 



E c 



2Cx 



(Qo + Nef 



(5) 



This is minimized when N is the integer closest to — Qo / e. 
Throughout this paper, we will assume that the Coulomb 
energy is much larger than the level spacing so that only 
the two lowest-energy values for N, namely No and N\ — 
Nq + 1, are permitted during the process of current flow. 
This assumption allows us to take the electrostatic energy 
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to be proportional to N: since 

is a constant, Eq for Nq or N\ electrons can be rewritten, 
to within a constant, as 



(a) 



(6) 



Notice that Eq. (|^) explicitly includes the Coulomb en- 
ergy which forbids states not having No or N\ electrons, 
but this is implicit in Eq. (^). The condition N = Nq 
or iVi has therefore to be assumed explicitly when using 
Eq. (§). 



2. Kinetic energy 

The kinetic energy of the electrons in the island can be 
written 



E K = '£efn i 



(7) 
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FIG. 2: Energy diagrams for the single electron transistor. 
The island is represented by a set of discrete energy levels 
and the leads by continua of levels. Filled dots in the island 
stand for electrons present in an iVo-electron ground state. 
The empty dot is an extra electron which tunnels onto the 
island to give an TVi-electron state. The transition marked 
with a solid arrow is the one which determines the initial 
threshold for starting current flow. The transitions marked 
by dotted arrows then also contribute to the total current, (a) 
When the Fermi energy of the right lead is swept past the first 
level available for tunneling at energy e^, current can tunnel 
through this level, (b) For a slightly lower gate voltage and 
higher bias voltage, two levels contribute to tunneling even at 
the initial onset of current flow. 



where ef is the energy, relative to the Fermi level, of 
spin-degenerate single-electron quantum state i, and m 
is the occupancy of this level (either 0, 1, or 2). 

Since N — J2i n i' tne sum °f the electrostatic and 
kinetic energies is just 



Eck = 



where is defined by 



K 



e 



Qq h o — e 



(8) 



(9) 



Writing the effective energy of the single-electron states 
in this way allows a simple accounting of the average 
Coulomb energy in the calculations. 

In the absence of variations in electron-electron inter- 
actions between electrons in different energy levels, the 
energy of the island is just Eck- With our conventions, 
the threshold voltages required for the onset of a tun- 
neling process can be pictured with simple energy dia- 
grams, as illustrated in Fig. g. For example, at T = 0, 
electrons can tunnel from lead k into the island if the is- 
land is a A^o-electron state and Fermi energy E£ of lead where td is the energy 
k is above the energy e; of a non-fully occupied level. 
In the same way, electrons can tunnel out of the island 
into lead k if the island is in a iVi-electron state and E£ 
is below the energy of a non-empty state. The onset of 
the current is associated with the first level available for 
tunneling, i.e. the lowest-energy non-full level in the Nq- 
electrons ground state or the highest-energy non-empty 
level in the ./Vi-electrons ground state. As V is ramped 



important to note that as V is increased, more than one 
spin-degenerate quantum level can contribute to tunnel- 
ing even at the initial onset of current flow. One example 
of this case is illustrated by Fig. ||(b) . The first allowed 
tunneling transition is for an electron to enter the level 
with energy from the right electrode. However, after 
this electron has tunneled in to give a total of Ni elec- 
trons on the island, transitions to the left electrode can 
occur either from the state with energy or from the 
lower-energy occupied state depicted in Fig. ||(b). If an 
electron tunnels out of the lower-energy state, subsequent 
tunneling transitions from the right electrode can involve 
either quantum level. Therefore, calculations of current 
for this situation must include tunneling processes occur- 
ring via both levels. 

It is possible to have current flow at vanishing V if the 
Fermi energy of both leads is aligned with the first level 
available for tunneling. The gate voltage that realizes 
this condition is called degeneracy point and is defined 
by 



V 



9i 

eC, 



(10) 



of this particular level. 



3. Variations in electron- electron interactions 

In the presence of variations in electron-electron inter- 
actions between electrons in different energy levelsjj the 
energy of the island has the extra term 



for a fixed value of V g , the Fermi energy in a lead can 
sweep past the energy required to initiate tunneling via 
an eigenstate, producing a stepwise change in current. 
The voltage-position, width, and current-amplitude of 
this step are the quantities that we will analyze. It is 



Ej = J({ni}) 



(11) 



Eq. (|l|) can be interpreted as an expansion of the energy 
of the system around the ground state: the second term 
is the part of E that is linear in {n^}, the first term is the 
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mean-field contribution of the quadratic part and J({ni}) 
is defined to be the rest. The net effect of the J({rii}) 
term is to produce shifts in the energy thresholds for 
tunneling that depend on the actual state of the particle. 
For instance, the effective energy level e[ for adding an 
election to level i starting with the A^-electron state {rij} 



Yi 



+ J({n j + S ij })-J({n j }). 



(12) 



Notice that this is only defined if rij < 2. In the same 
way, the energy of a non-empty energy level in a N\ state 
can be defined as minus the energy required to remove 
an electron from that level. 



B. Steady-state occupation probabilities 

Because of the influence of the Coulomb charging en- 
ergy, even in the simplest cases that we will consider the 
occupation probability for a given many-body state \a) — 
{n-i} of the particle cannot be factorized as the product 
of occupancy probabilities for each single-electron level. 
Therefore we have to solve the full rate-equation prob- 
lem where the occupation probability of each many-body 
state is treated as an independent variable. 

The evolution, of the occupation probability of state 
| a) is given byoQ 

-f = "£(r^ a Pp-T a ^p a ) (13) 

p 

where T a -^$ is the transition rate from state \a) to state 




FIG. 3: Energy diagram with one level available for tunneling. 



is the Fermi function corresponding to the temperature 
in the leads and 7- and 7J are the bare tunneling rates 
between level i and each of the leads. Here e[ is the 
energy needed to add an electron to state \a) in level i. 
It includes the contribution of the interaction term. 

The steady-state occupation probabilities can be found 
by iterating Eq. ( |l4| ) with a discrete timestep dt to find 
the probabilities for which dP /dt = 0. This is equivalent 
to finding the eigenvector P of T associated with the 
eigenvalue zero. 



Current 



Once the occupation probabilities for each state |a) are 
determined at given values of V and V g , then the current 
can be calculated either through the right tunnel barrier 
or through the left barrier. In the steady state these two 
cyj^ents are equal. The current through the left barrier 



(18) 



This can be written in matrix form as 
dP 
~dt 

with the following coefficients for the matrix T: 



= r-p 



1 a/3 — 

r — 

A aa — 



P' 



if 



P^a 



(14) 

(15a) 
(15b) 



We do not consider cotunncling or internal relaxation 
in the particle. Therefore, the only states that are cou- 
pled together are states that have the same occupancy 
for all the levels, except one electron difference in one 
level. Let's assume that states \a) and \f3) differ only by 
having one extra electron in level i. Then 



where 



r a ^ = 71/(^-^x2-0 

+ 7[/(^-Sf)(2-n i ) 
T^ ot = 1 \(l-f{e l i -E())n i 

f{x) = 1/(1 + exp(x/k B T)) 



(16a) 
(16b) 

(17) 



where T a ^g is the contribution of the left lead to r a ^>p, 
multiplied by + 1 or — 1 depending on whether the a —> f3 
transition gives a positive or negative contribution to the 
current. 

In order to get a feeling of the physics that will come 
out of this rate-equation model, in the rest of the paper 
we will consider selected examples that are simple enough 
to be solved by hand, yet have the basic ingredients of 
the complete problem. 



III. ONE SPIN-DEGENERATE LEVEL 
ACCESSIBLE 

A. General formula 

Consider the situation represented in Fig. ^ where only 
one spin-degenerate energy level, with energy e\ , is acces- 
sible for tunneling and (on account of the large Coulomb 
energy) it cajx.be occupied by either zero or one electron, 
but not two a If we call 



fr = f(ei 
fi = f(ei 



eh 



(19a) 
(19b) 
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FIG. 4: Current profiles as a function of the bias voltage 
for the case of a single spin-degenerate level accessible for 
tunneling, for three different gate voltages. We assume CV — 
Ci and ji — 4^y r . The bias voltage is plotted in units of 
fcsT/|e|. The current is in units of |e|-yo where 70 = 7;7r/(7; + 
7,.). The reduced gate voltage v g = |e|C 9 (V 9 - V°)/C s k B T 
is 0, -3 or -6. 



and N — or 1 the state with N electrons, the transition 
rates are 

IW = 2 lr f r + 2 7i / ; (20a) 

ri_ >0 =7r(l-/r)+7j(l-/l) (20b) 

for the tunneling-in and tunneling-out transitions. Then, 
the occupation probabilities are 



Pi = 
Pa = 



2 7r / r + 2 7i / ( 



7r(l + /r)+7/(l + /0 
7r(l-/r)+Ti(l-/l) 



(21a) 
(21b) 



IW+r^o 7r(l + /r)+7l(l + /l) 

and the current through the left lead in the steady state 



/=|e|(7 / (l-/ / )P 1 -2 7/ / / Po) 

7r7i(/r - //) 



2|e| 



7r(l + /r)+7*(l + /0' 



(22) 



This expression differs from an approximate form used in 
Ref. ^ to analyze tunneling data. 

We can plot the current as a function of the applied 
voltages by replacing fk by their definitions in Eqs. jig| ) 
and by the expressions in Eqs. (|J). Fig. |J shows the 
current steps as a function of the bias voltage when the 
gate voltage is first equal to the degeneracy point, then 
is tuned away from it. 



B. High bias limit 

If the level spacing is very large compared to fegT, 
there is an interesting regime in which V is substantially 
bigger than fcgT/|e| yet only one level is involved in the 
current transport. The limiting current in this case is 
bias-independent and can be obtained from Eq. ( p2] ) by 



FIG. 5: Energy diagram with one level available for tunneling 
and V g < V§. Since Ef is substantially below ei, electrons 
can tunnel into the island only from the right lead. 



setting f r = 1 and fi=0 (positive bias) or // = 1 and 
f r = (negative bias). For these two cases we have 
respectively:!^ 



2\e\ 



-2e 



TrTj 

27r + ll 
Irll 



7r + 27; 



(23a) 
(23b) 



These expressions give different heights for the positive 
and negative current steps. Measuring these heights can 
therefore allow an experimental determination of both 
7 r and 7;. Note that this is in contrast with the case in 
which tunneling occurs through a single level that is not 
spin degenerate. In that case 



it = ±M 



7r7j 
7r +7 



(24) 



for both bias directions ,0 so that 7 r and 7; can not be 
determined separately. 

In the limit of barriers with very different tunneling 
rates (which can be experimentally relevant if the barrier 
thickness is not well controlled), the current depends only 
on the smaller 7. For example, if 7; 3> "/ r , then I + = 
2|e|7 r and L_ = — |e|7 r . The factor of 2 in I+/I- arises 
from the difference in the number of spin states accessible 
for tunneling for the rate-limiting transition across the 
right barrier. 



C. Position and width of the current step 



Next we consider the case depicted in Fig. |^, in which 
V g is adjusted away from the degeneracy point so that 
at the threshold V for tunneling only the effective Fermi 
energy in the right electrode is close to e\ , while the Fermi 
energy of the left electrode is at a much lower energy. 
That is, we will assume /; = 0. Using this assumption, 
after some algebra Eq. (p3) becomes 



/ = /+/ ei-ST -fc B Tln 



27r + 7 
7r + 7 



(25) 



Even though both spin-states of the quantum level con- 
tribute to tunneling, we can see in this expression that 
the current step has the shape of a simple Fermi func- 
tion whose width is given by the electron temperature of 
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FIG. 6: (a) Current step and (b) conductance peak at posi- 
tive bias and negative gate voltage for three different temper- 
atures. We assume C r = C;, 7i = 50 MHztl and j r ~> 7;. The 
peak occurs at Vb = 2C g (V g — V^°)/Cs at zero temperature 
and shifts from this position by an amount 2fesTln2/|e| at 
non-zero temperature. 



the leads. However, at non-zero temperature, the cen- 
ter of the step is shifted relative to its position at zero 
temperature. The shift is proportional to the temper- 
ature, vanishes if 7; 3> "f r , and has a maximum value 
of fcsTm2 when j r ^ 7/. Fig. ^ shows the shape of 
the conductance peak dl/dV in the latter limit for three 
different temperatures. 

There is a simple intuitive explanation of the shift in 
the limit j r 3> 7/ . The current threshold at zero temper- 
ature is given by = e\. At non-zero temperatures, 
when = ei, the Fermi occupancy probability is 1/2 
for states in the right lead with energy ei. In this case 
the transition rates are dominated by electrons tunneling 
back and forth from the right lead (since j r 3> 7;): 



r - 



r 



(27,) 



1^0 



(26a) 
(26b) 



D. Zeeman splitting of the energy level 

In the presence of an applied magnetic field, the two 
spin states associated with a given orbital level are no 
longer degenerate, but split to give the energies = 
£1 ± gfj.B^oH/2. If we call these states + and — , and 



E£), then the transition rates are 



lift 



T ^± = Irft 

r ± ^ = 7l .(l-/ r ± )+7;(l-/ J ± )- 



(27a) 
(27b) 



Notice the absence of the factors 2 that were in Eqs. ( |2(j| ) 
due to the spin degeneracy. The occupation probabilities 
are 



1 



1 



... I 



r ^± p 

-"a 



and the current through the left lead is 



(28a) 
(28b) 



i = ie| 7 /((i - f+)p+ + (1 - mp- - (ft + fn p o) 

(29) 

Figure [?] shows the effect of the magnetic field on the 
conductance peak at positive bias for a gate voltage be- 
low the degeneracy point (i.e. the case ff = 0). The 
peak splits into two subpeaks of different weight. This 
asymmetry can be understood by noticing that the first 
subpeak carries a current given by Eq. ( |24| ) and the two 
peaks together give a total current given by Eq. (|2^). 
Then the fraction of the total current carried by the first 
subpeak is just 



If 



i£ = 2 7r + 7; 

7+ 2 7r + 2 7; ' 
this ratio is one and 



(30) 



7 r ^—.7;, tins ratio is one ana the second peak 
vanishes. E£l On the other hand, if 7; 3> j r , the peak splits 
into two subpeaks carrying the same current. 



Here the factor 1/2 comes from the Fermi occupancy of 
the lead and the factor 2 in r ^i is from the spin de- 
generacy. This factor is only present for To^i because 
electrons tunneling into the island see two empty states, 
while an electron tunneling out comes from a given spin 
state. It follows that when Ef = e\ the probability that 
the island is in the 1-electron state is exactly two-thirds. 
Then the rate at which electrons tunnel to the left lead 
(the rate-limiting process determining the total current) 
is two-thirds of the maximum value. This can be seen 
directly in Fig. ||(a), or in Fig. ^(b) by the fact that two- 
thirds of the current (area under the peaks) lies left of 
V — Vq = 0. This T-dependent shift in the apparent 
resonance position has been observed by Deshpande et 
alE 



IV. TWO LEVELS ACCESSIBLE 

Next consider the situation pictured in Fig. |^ where 
two spin-degenerate levels are accessible for tunneling 
and the number of electrons in these levels is N — 2 or 3. 
Because of Coulomb blockade, no current flow is possible 
until an electron can tunnel from the right electrode to 
state 2; however after this happens both states 1 and 2 
can contribute to the current even at the initial current 
onset. Let (n\, n 2 ) be the state with ri\ electrons in level 
1 and n-i electrons in level 2 (n\ +ri 2 = TV), let P(ni,n 2 ) 
be the probability of state (m,n 2 ), and let 7^' be the bare 
tunneling rate of level i across barrier k. We will special- 
ize immediately to the interesting case of positive bias 
(as pictured in Fig. |^) with the right barrier substantially 
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FIG. 7: Splitting of a conductance peak in a magnetic field. 



We assume /, 



0, C r = Ci, and y r 



1i 



7- 



At zero 



temperature and zero field the peak occurs at Vb = 2C g (V g — 
V°)/Cs- The reduced field h = g^ B ^oH/(2k B T) is 0, 3 or 6. 




FIG. 8: Energy diagram for a case with two levels available 
for tunneling. 



thicker than the left barrier, so 7' ^> 7 J for i,j £ {1,2}. 
(Note that this is opposite to the inequality considered in 
Fig. ||.) To simplify further, we will also look only at the 
current onset at positive bias for a large negative gate 
voltage, i.e. we will assume f(e\ — Ef) = f(e 2 — Ef) = 
and f(ei-E^) — 1. For this case, f(e 2 — E^) will simply 
be called /. These conditions correspond to line III in 
the data of Ref. ||. 

Figure || shows the available transitions together with 
the corresponding transition rates. Since <C 7? for 
i = 1,2, the terms having a factor 1 — / can be neglected. 



A. Rate equation 

The rate equation in this case has to describe eight pos- 
sible transitions between five different states. It is there- 
fore convenient to use the matrix notation of Eq. (pd|), 
which gives 



dt 
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p(i,r x 

P(0,2 
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P(2,l) 



(31a) 
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-2/ 72 r 








li + (1 - fH 










-il - hi 





H 


H + 2(1 


- f)l2 








-27l 









2/7 2 r 


Yi 





-l l 2 - (1 - f)l2 ~ H 







, 


hi 


27[ 





~H - 2(1 - 


fH - -A 



(31b) 



This matrix has the structure 



B. Current 



r = 



r ^ uu 




r cu 


Tec 



(32) 



where T uu and r cc are diagonal blocks associated respec- 
tively with the 7V -electron (uncharged) and A^-electron 
(charged) states. The cross-diagonal blocks are associ- 
ated with the tunncling-out (r uc ) and tunncling-in (r cu ) 
events. This structure is preserved whatever number of 
levels are available for tunneling. 

In the steady state, the solutions for the occupation 
probabilities are as follows: 



P(2,1)<1 P(1,2)<1 



P(2,0) = - P(l,l) 



AfK 
S 



P(0,2) 



f 2 K 2 
S 



(33a) 
(33b) 



where K = 2fa and S = 1 + ifK + f 2 K 2 . 

T2T1 



Since we can neglect the tunneling-out transitions 
through the right barrier, we can calculate the current as 
the sum of the contributions of the tunneling-in events 
through this barrier: 



2fY 2 P(2, 0) + (/ 72 r + Y)P(l, 1) + 2 7 [P(0, 2) 
. (4Y 2 K + 2 7 [X 2 )/ 2 + (2Y 2 + 4 7 D/ 



1 + AfK + f 2 K 2 



(34) 



In Fig. |l0| we compare this expression to the current we 
would have in the presence of infinitely fast relaxation in 
the island (state (1, 1) relaxing instantaneously to (2, 0)). 
In such a case electrons can only tunnel into the higher 
energy level in the island. Since the tunneling-in of elec- 
trons is the rate-limiting process, this situation is equiv- 
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FIG. 9: Available transitions for the situation described in Fig. ^. 




|e|(V - V )/k B T 



FIG. 10: Shift of the current step by non-equilibrium in the 2- 
levels- accessible case. We assume 7} = 72 = 7; and 7i = 72 = 
j r , with 7; 3> 7r- The step occurs at Vo = 2C g (V g — Vg)/Cs 
at zero temperature. The "equilibrium" curve assumes in- 
finitely fast relaxation in the island. The "non equilibrium" 
curve assumes no relaxation. 



alent to the case of Eq. ( |22| ) when only one level is acces- 
sible for tunneling, and the current would just be 
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FIG. 11: Dependence of (a) the current step and (b) the con- 
ductance peak on the temperature in the 2-levels-accessible 
case in the presence of non-equilibrium. We assume 71 = 
72 = ji and 71 = 72 = 7r, with 7; ^> j r and no relaxation in 
the island. 



^equilibrium = 2\e\"f2 f ( e 2 ~ 



(35) 



The main effect of non-equilibrium states as illustrated 
in Fig. [n] is therefore to shift the current step to lower 
voltage. Although not exactly a Fermi function, the 
shape of the step described by Eq. ( |34| ) is very close to 
a Fermi function, shifted by — 1.79fcsT and widened by 
8.5%. The shift can be understood as follows: When 
Ef = £2, electrons tunneling to the upper level come 
from half-full states in the right lead. If the island is in a 
non-equilibrium state ((1, 1) or (0, 2)), electrons can also 
tunnel to the lower level. Since these electrons come from 
full states in the lead, the current at Ef = e 2 is higher 
when these states are allowed, hence the shift. 

The temperature dependence of the current step and 
the conductance peak in this two-level-accessible case 
with 7; ^ 7,. is displayed in Fig. [ll]. Although the In- 
dependent shift looks very similar to the result for one 
level displayed in Fig. ||, the shift in Fig. [ll] is of a differ- 
ent nature since it originates from non-equilibrium states. 
For the one-level-accessible case, there was no shift for 
positive bias with 7; 3> "f r . If we look at the opposite 
limit with two levels (positive bias 7; <C 7 r ), the rate 
equation will be dominated by electrons tunneling back 
and forth between the right lead and the second level in 
the island. This situation is very similar to the one-level 
case and gives the current: 

I = \e\(2j[ + l2 )f(e 2 - E F r - fc B Tln2), (36) 



where the shift by fc^T In 2 is explained by the same argu- 
ment as in the one-level case. The additional level there- 
fore does not produce an additional shift when 7/ <C 7r- 

If the voltages are tuned so that more than two levels 
are made available for tunneling-out transitions (by low- 
ering Ef), or if the tunnel couplings to state 1 are greater 
than to state 2, then the shifting of the resonance away 
from the T — position in the 7; 3> tv case will be en- 
hanced beyond what is shown in Fig. nil This shift will 
however remain proportional to fceT. 

We have also considered the case when Ef is very high, 
so that many levels are accessible for an electron to tunnel 
into the island across the higher-resistance tunnel bar- 
rier, while Ef remains fixed slightly below e-y. In such 
a situation the total tunneling-in transition rate will be 
proportional to the number of levels available for tunnel- 
ing in, and this rate can eventually become greater than 
the tunneling-out rate which will be roughly constant. In 
this case, tunneling though the left lead will eventually 
become the bottleneck process even if 7/ » j r) which 
allows one to estimate an average tunneling rate through 
the lower-resistance barrier even in the case of very asym- 
metric barriers. 
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FIG. 12: Current steps for different interaction strengths. We 
assume C\ = C r , j{ = 72 = 7i, and 7! = J2 = Jr, with 7; > 
7 r . The "equilibrium" curve assumes infinitely fast relaxation 
in the island. The other curves assume no relaxation and 
S/ksT ranging (a) from to -20 and (b) to 6. 



V. TWO LEVELS ACCESSIBLE WITH 
VARIATIONS IN THE INTERACTIONS 

In the presence of variations in electron-electron inter- 
actions, the energy thresholds for tunneling are different 
depending on whether the island is initially in a ground 
state or in an excited state. For example, in the case de- 
scribed in the previous section, this effect can make the 
energy required for the (1,1) — > (1,2) transition differ- 
ent than the (2,0) — > (2, 1) transition. We can account 
for such variations by assigning a different energy to the 
upper level in the presence or absence of an excitation 
in the island. Namely, the energy of the upper level will 
be £2 for the (2,0) — > (2,1) transition and e 2 = €2 + 8 
for the (1,1) — > (1,2) transition. Here 8 is a measure 
of the strength of the variations. In order to generalize 
the previous notation, we will call / = f(e2 — Ef) and 

f' = f&-E?). 

The possible transitions are still described by Fig. ^ 
and the corresponding rate equations are the same as 
Eqs. <M) but with 



r 



(1,1)^(1,2) 



r ( i, 2 )^(i,i) = 2 72 + 2(1 - /') 72 r 
which gives the current 

/= , (4 72 r ^ + 2 7 r^ 2 )// ; + (2 72 r + 4 7 D/ 
11 1 + ifK + ff'K* 



(37a) 
(37b) 



(38) 



Fig. [12](a) shows the current step for the case that the en- 
ergy required for the tunneling transition is decreased by 
non-equilibrium (negative 8) for various values of 8/ksT 



ranging from to -20, and Fig. |12|(b) shows I-V curves 
when the non-equilibrium effect increases the tunneling 
energy. These plots were made for 7]" = 72 = 7r , 7/ 3> 7r , 
and K = 1. For negative 8, the effect of the variation 
in electron-electron interactions is to produce an addi- 
tional shift in the voltage-position of the current step, 
on top of the shift already described due to the non- 
equilibrium states. This additional shift is proportional 
to 1 8 1 if 1 5 1 <C ksT and becomes a constant on the order 
of k B T if \S\ 3> k B T. A shift of this sort has been ob- 
served in Fig. 3(b) of Ref. |8|. For positive 8, the effect of 
non-equilibrium is to produce an extra step in the I-V 
curve at voltages larger than the position of the (5 = 
current step. 

As V is increased so that more than two levels be- 
come energetically-accessible for tunneling, the ensemble 
of possible non-equilibrium excitations grows combina- 
torically, and each combination of excitations can pro- 
duce a different shift for the tunneling resonance ener- 
gies. Interactions which depend on the spin state of the 
island (neglected thus far) can produce further compli- 
cations. The non-equilibrium excitations can produce a 
variety of effects depending on the ratio 7 ;/7 r and on 
the magnitude of variations in electron-electron interac- 
tions. When the interaction-induced shifts are compa- 
rable to fcflT, they have been observed to produce an. 
effective broadening of the observed conductance peaks .El 
For larger interactions, shifts due to non-gauilibrium ex- 
citations have been resolved individually.LrE3 



VI. CONCLUSIONS 

We have solved the rate equations describing electron 
tunneling via discrete quantum states on a nanoscale is- 
land, for selected simple cases, under the assumption that 
rate for internal relaxation of excited electronic states is 
slower than the electron tunneling rate. Even the sim- 
plest case of tunneling via a single spin-degenerate energy 
level has some initially-surprising features. The magni- 
tude of the maximum tunneling current can depend on 
the sign of the applied bias V, and the voltage-position 
of the resonance is temperature-dependent. When two 
spin-degenerate quantum levels are accessible for tun- 
neling, the behavior is even richer because of the influ- 
ence of non-equilibrium excitations on the island. The 
voltage-position of the resonance can undergo strong 
temperature-dependent shifts even in regimes (e.g., pos- 
itive bias and 7 / 3> j r noted above) where the one-level 
resonance positions do not depend on temperature. Un- 
derstanding the variations in the strength of electron- 
electron interactions is critical in the non-equilibrium 
regime with two or more levels accessible. Such varia- 
tions can produce additional shifts of resonance curves 
on top of the shifts noted previously, and they can also 
introduce extra steps into the current- voltage curves. 

The methods we have described for determining tun- 
neling currents are applicable to more than two levels, 
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but the analytic expressions become sufficiently compli- 
cated to be of limited usefulness. We have verified nu- 
merically that the results for additional levels are qualita- 
tively similar to the two-level case. The computer codes 
we have used for calc ulating the general cases are avail- 
able elec tronically at \ittp : //www, ccmr . Cornell . edu/ 



"ra.ph/. These are useful, for instance, in extracting 
the rate-limiting bare tunneling rates from experimental 
data in which stepwise increases in current are measured 
as V and V„ are adjusted, so that the number of states 



accessible for tunneling increases one-by-one.l 
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